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Abstract 

<N : 

We prove a new soft pion theorem for the near threshold pion production by a hard 
electromagnetic probe. This theorem relates various near threshold pion production 
Ph ■ amplitudes to the nucleon distribution amplitudes. The new soft pion theorem is 

t^. | in a good agreement with the SLAC data for the structure function F^iW, Q 2 ) for 

W 2 < 1.4 GeV 2 and 7 < Q 2 < 30.7 GeV 2 . 
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Introduction 

(N ; 

i [ The amplitudes of the pion production near the threshold by the electromagnetic probe 
O . 

g ; 7* + N -> 7T + N' (1) 

at a not too large virtuality (Q 2 ) of the photon Q 2 <C A 3 /m 7T (A ~ 1 GeV is a typical 
hadronic scale) can be expressed in terms of various nucleon form factors with the help of 
the soft pion theorem (SPT) [[I], |^, [|. For virtualities Q 2 ~ A 3 /m n and larger this SPT does 
jg: not work 0,|, | . 

In the present paper we derive a new "hard-soft" pion theorem (hSPT) for the reaction 
([]]) in the near threshold region and for large virtuality of the photon Q 2 ^> A 2 . Our main 
tool is the QCD factorization theorem for exclusive processes {|, || (for a recent review and 
comprehensive list of references see e.g. ||). It allows us to express the pion production 
amplitude at large virtuality in terms of the distribution amplitudes (DAs) of the nucleon 
and of the low-mass ttN system. These non-perturbative objects correspond to the lowest 
three quark (3q) Fock component of the nucleon and ttN systems. We derive a hSPT to 
relate the corresponding distribution amplitude of the 7rN system (we shall call it as tcN DA) 
to the nucleon distribution amplitude. This hSPT is valid for the limit when the mass of the 
ttN system (denoted as W) is close to the threshold W t h = M N + m^, i.e. W — W^^m^. 
The derivation of the similar theorem for DA of the two pion system near threshold can be 
found in 0,| . 

The physical picture of the near threshold production of pion by a hard electromagnetic 
probe is as follows. At large Q 2 the emission of the soft pion from the initial state contributes 
only to large invariant masses W . The emission from the hard interaction part is a higher 
twist in Q 2 . Hence the emission occurs solely in the final state when a small 3q system 
produced in the hard interaction expands to a large enough configuration. At this point 
we are dealing with a soft pion emission and can apply corresponding near threshold chiral 
theory relations. 
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Threshold theorem for hard electromagnetic probe 

We consider the hard reaction ([I]) in the near threshold region, i.e. when W is close to the 
threshold of the reaction W t h = + m n and at large virtuality of the photon Q ;§> A. 
Hence, xsj — »• [1 + m n (2MN + rn^) / Q 2 } 1 in the discussed limit. In the near threshold region 
the cross section of the reaction (jl]) can be expressed in terms of the N —* ttN transition 
form factors at large momentum transfer Q 2 : 



Mf/ a = (T*N f \JZ*-{G)\N f ) (2) 

Here /, /' = p,n are flavours of the initial and final nucleons, a is the flavour index of the 
emitted pion. 

In this section we shall be interested in the matrix element (|2|) at the threshold of the 
reaction (|l|), i.e. at W = W t h- The matrix element (0) at the threshold for the transverse 
photon has the formQ: 



Ml{ a = A(j*Nf - N f , + vr a ) w(p')7^76«(p) • (3) 

It follows from the QCD factorization theorem {|, |5j that the transition matrix element 
A(j*N — > ttN') at large Q 2 can be written as (up to the power suppressed terms) 

A ^*N - ttN') = J dxdy<S>l N ,(x)T(x,y,Q 2 )® N (y) , (4) 

where T(x, y, Q 2 ) is the hard part of the process computable in perturbative QCD. The 
functions ^^(y), Q n Ni(x) are distribution amplitudes (light cone wave functions) of the 
nucleon and of the low- mass ttN' system. The DA of the nucleon <$>n{u) a ls° enters the QCD 
description of the nucleon form factor and it is a subject of intensive studies. The distribution 
amplitude of the ttN system is a straightforward generalization of the baryon DAs. Many 
general properties of these objects are similar to those of the two-pion distribution amplitude 
which were extensively studied in Ref. ||. They will be discussed elsewhere. 

We focus here on the soft pion theorem for irN DAs. Using the general soft pion theorem 
(see e.g. 0) we can write: 

(7r a (k)N f (p,S)\O\0) = UN f (p,S)\ [Q a 5 ,0] |0> (5) 

Jtt 

+ Af X. k) 5>(P> S)HT] f ,u(p, S')(N f ,(p, S')\O\0) . 

Here f n = 93 MeV is the pion decay constant, ~ 1.25 is the axial charge of the nucleon. 
The operator O is the trilocal quark operator of twist-3 in which quark fields are separated 
by a light cone distance, in the complete analogy with the definition of the baryon DA, see 
e.g. (Til. The operator Q§ is the operator of axial charge, so that the calculation of the 



commutator [Qf, O] reduces to the chiral rotation of operator O. The second term in the rhs 
of Eq. (HI) corresponds to the chiral singularity due to the nucleon pole in the graph shown 
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Figure 1: Nucleon pole contribution to the soft pion theorem for generalized ttN distribution 
amplitudes. 

in Fig. 1. The contribution of this diagram is strongly suppressed for W — W t h <C ra % but 
for W — W t h ~ Tfi-K it becomes significant, see Eqs. (|2^j25|) . 

Let us start from the calculation of the first (commutator) term in Eq. (|5|). Since the 
chiral rotation of the trilocal quark operator O does not change its twist Eq. fl5|) allows us 
to express generalized ttN DAs at the threshold in terms of nucleon DAs. 

We write the nucleon DA in terms of functions (fis( x ) anci 4>a(x) which are symmetric 
and antisymmetric respectively under x\ <-> x% (1 and 3 are quarks with parallel helicities) 
0) Hi- For ^ e P r oton we have 



+ ^^NMt- rf T M i M t>- (6) 

The distribution amplitude for neutron can be obtained from the above expression by the 
replacement u <-> d. 

Applying the general soft pion theorem (|5]) we express the distribution amplitudes of 
various ttN systems at the threshold in terms of the nucleon DAs <fis{ x ) an d 4>a (%)'■ 

\P 1 7r °) = t S ^} I 6m t^i m t + «t«|d T + d^u\u^) 



(Pa(x) 



u 1 u l d 1 - ^t m i m t) > ( 7 ) 



|n t 7r+ ) = *tjl±!L- |2w-[-dj_u-f — 3 u^u^ — ^ d-fUjii-]-) 



2/. 



|w T W;d T - d T W> • ( 8 ) 



The DAs of the neutral ttN systems can be obtained by the trivial replacement u <-> d. 

Now we can compute the threshold amplitudes A{^*N — > ttN') at large Q 2 combining 
the factorization theorem (^) with the expressions for nN DAs The technique of 

calculations of the hard part T(x,y,Q 2 ) is standard and can be found, e.g. in refs. ||, [TO]. 

* Longitudinally polarized photon gives a power suppressed contribution for Q 2 3> A 2 . 
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We give here only the final results using notations of refs. j|, |TT|, [E^. For the transitions 
p — > 7riV we obtain: 



th 



(167ra 8 )- 
9/. 



| [dx dy] {^(23Tx - 8T 2 ) <j> s (x)Mv) (9) 



1 



v/3 



<Psi.x)(p A {.y) + <Ps(v)<Pa(x) 
- ^T 1 + 2T 2 )0 A ( 2/ )0 A (x)|, 



Q 4 A(7*p -> 7T+n) 



th 



y [dx dy] {^(lm + 4T 2 ) (10) 

(j)s{x)^A{y) + <Ps(3j)<Pa{.x) 



2 T 

7! Tl 



T 1 + 2r 2 U A (y)0 A (x) . 



Here [dx] = dx\dx 2 dx 3 5(1 — x\ — x 2 — 23) and coefficient functions T% and T 2 have the form 



1 



x 3 (l - x^y^l - yi) 2 x 2 (l - x 1 ) 2 y 2 {l - yi) 2 
1 
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x 2 x 3 (l - x 3 )y 2 y 3 (l - Vi) 
1 

£1X3(1 - 2:1)3/12/3(1 - 2/3) 



(11) 



(12) 



Expressions for various nucleon form factors in the same notations can be found e.g. in 
refs. § 0, 0- 

For the transitions n — > ttN we obtain: 



Q A A(-f*n -> rr°n) 



th 



I [dx dy] {^(t\ - r 2 ) M*)My) (is) 



(167ra a ) ; 
9/. 

+ -7=^1 <Ps{x)<pA(y) + <t>s(y)(f>A(x) 
+ ±(r 1 -T 2 yA(y)Mx)}, 



Q 4 A( 7 *n^n-p)\ 



th 



(16vro; s ) : 



j[dx dy] j-^Tx-Ta) 



MtfMv) (14) 



+ 7S T ' 



<f>s{x)(f>A{y) + <f>s{y)<f>A{x) 



Tx - T 2 )0 A (2/)0 j4 (a;) 



Using Eq. (Q) one can express the near threshold structure functions Ff ,n (W, Q 2 ) and 
various differential cross sections for particular ttN channels in terms of the nucleon DAs. 
The same DAs appear in the QCD factorization theorem for the nucleon form factors at 
large Q 2 . We find that in the case of the symmetric form of the nucleon DA 

<ps (x) is arbitrary, 0^(2;)= 0, (15) 

one can describe the near threshold pion production directly in terms of the nucleon form 
factors without specifying the nucleon DA, see below. Certainly the nucleon DA can have a 
nonzero asymmetric component <pA Actually our general Eq. ([]) can be applied to any 
specific model of the nucleon DA (see e.g. [[C|) to compare the model predictions with the 
experiment. In the case of symmetric nucleon DA the amplitude of the process (|l|) can be 
expressed in terms of nucleon magnetic form factors {Gmn{Q 2 )) as follows: 

A(7*p-7r°p)k = -j (J G *«p ~ t G ^ > ( 16 ) 
A( 7 *p -> 7r+n)| th = -i- (J G Mp + \ G Mn ) , (17) 



A(j*n - ir°n)\ th = (18) 

6/vr 

Atfn^TT-p)^ = (19) 

Since these results are based on the symmetric form ([T^) for the nucleon distribution ampli- 
tude, deviations from these equations would allow to probe directly the asymmetric part of 
the nucleon distribution amplitude and check the validity of the leading twist description of 
the nucleon form factors. 

At asymptotically large Q 2 the DAs of the nucleon are evolved to their asymptotic values: 

C(x) = N x lX2 x 3 , <^ s (x) = . (20) 

As the form (up to overall normalization) of the asymptotic nucleon DAs is fixed uniquely 
we can make predictions for the ratios of various amplitudes at the threshold at Q 2 — > oo 

Zz&t = o(flt). Z^Sn = -« + °{m) 

are qualitatively different from predictions for the case — > JT]] limg2^ ^^ZTJ+n)^ = 
O (^jjf^j F° r the ratio of neutron to proton structure functions we obtain from new soft pion 

theorems the following asymptotic result: liniQ2_> 0O)W '->W th ^ w '^ — |^ + Ol j^) , which 

is much larger than the perturbative QCD scaling limit expectation of 3/7 |L4j]. Although it 
is worth noting that the presently experimentally accessible values of Q 2 are not very large. 
Therefore one may expect considerable deviations from the above asymptotic results. 
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Structure functions: comparison with data 

The data on the structure function i*-f(W, Q 2 ) in the near threshold region were obtained in 
1994 in the SLAC experiment E136 [l5j for a wide range of Q 2 = 7-j-30.7 GeV 2 . Here we make 
the first comparison of the data with the hSPT predictions. (Previous analyses addressed 
the scaling features of these data and did not attempt to calculate F$ {W < 1.2GeV, Q 2 )-) 

To compute F^iW^Q 2 ) for W — W^^m^ we combine the strictly threshold amplitude 
@ with the contribution of the last term from Eq. ([5]) and obtain 



•£ \A(;-p ^ X)U (21) 

X=piz°, nir+ 

2^ /32/OAU/4 



3g 2 A G 2 Mp (Q 2 )(3 2 (W)W* (m « 
A{W 2 -M 2 N + ml) 2 V A 



where 



0m = ^«^L^(«^. (22) 

The first term in the rhs of eq. (]2l|) corresponds to strictly threshold amplitude (||). The 
second term takes into account the emission of the pion from the outgoing nucleon (Fig. 1) 
with the amplitude given by the second term in the rhs of eq. @. The latter contribution 
vanishes exactly at the threshold but gives a parametrically unsuppressed contribution for 
W — Wth ~ m 7r- Note that the second term in Eq. corresponds to ttN system in P- 
wave, therefore it can be separated from the first (S-wave) term by considering the angular 
distributions in ttN system. 

The data of E136 experiment Ref. |D| are consistent with the factorized form 

Fi(W,Q 2 )=F(Q 2 )G(W), (23) 

with F(Q 2 ) oc l/Q 6 for Q 2 > 8 GeV 2 which is exactly the scaling behaviour following from 
Eq. fl21|). We also found that Eq. (|2T|) provides a good description of the W-dependence of 
the E136 data for W < 1.2 GeV though the predicted jy-dependence is somewhat different 
from the G(W) oc W 2 — W 2 h fit of However, the resolution of E136 is not sufficient 



to distinguish between the two forms. The hSPT (|2i|) also predicts the scaling behaviour 
~ 1/Q 6 which is confirmed by the E136 data [jTjJ. 

Thus we reproduce several features of the data without using any specific nucleon wave 
functions. To make a first quantitative comparison of the soft pion theorem prediction for 
the absolute value of F$ we use the symmetric form for the nucleon distribution amplitude 
(|j~5|). Inserting the expressions ( |i~6|) and ( |TTD for the threshold amplitudes into Eq. (|21"D and 
using the analogous expression for the neutron structure function, we obtain the following 
hSPT for the structure functions in the case of the symmetric form ( [T5D for the nucleon DA 



+ 



WW [25 2 2 
[l2 Mp [Q } 
3g 2 A G 2 Mp (Q 2 )f3 2 (W)W* 
A(W 2 — M% + ml) 2 



+ o G 2 Mn (Q 2 ) 



(24) 



+ 



rrh 
A 



6 



F?(W,Q 2 



Q 2 (3(W) 



— G 2 
4 



Mn 



(Q 2 



+ 



Zg 2 A Gl In {Q 2 W{W)W± 
4{W 2 



(25) 



M 2 N + ml] 



+ 



A 



Now we can compare our results based on the hSPT with the E136 data ||15| . It was found 
in this experiment that in the interval W 2 < 1.4 GeV 2 the quantity J^ 4 dW 2 Q 6 F%(W, Q 2 ) 
is practically constant for Q 2 > 8 GeV 2 . The experimental value of the integral is in a good 
agreement with our result 



1.4 



dW 2 Q 6 F%{W, Q 2 



th 



0.10 ±0.02 GeV 8 (E136) 
0.11 ±0.02 GeV 8 (hSPT) 



(26) 



For the theoretical analysis we used as input the following values for the nucleon form factors 
Q 4 G Mp (Q 2 ) = 1-0 ± 0.1 GeV 4 obtained at Q 2 > 10 GeV 2 and Q 4 G Mn (Q 2 ) = -(0.5 ± 
0.1) GeV 4 extracted from Ref. |17| at Q 2 ~ 10 GeV 2 fll8 |. For the summary of the current 
information on various baryon form factors at large momentum transfer see e.g. [|19J. Note 
that the contribution of the P-wave term in Eq. (|24]) is relatively small (about 20%), so that 
the main part of the hSPT value (p6|) is due to the strictly threshold contribution (S-wave) 
in Eq. (|l|)- With the same set of parameters we predict J^ 4 dW 2 Q 6 F£(W, Q 2 ) = 0.05 ± 
0.02 GeV 8 . Note that for the nucleon DAs which fit G Mp (Q 2 ) , G Mn (Q 2 ) at Q 2 > 10 GeV 2 
(such DA significantly differ from the asymptotic one, see e.g. |ll], [L2|]) the F^/Frj ratio near 
threshold is much smaller than the asymptotic value of 52/37 which follows from Eq.(^) with 
the asymptotic distribution amplitude <ft(x) oc £1X2^3- Thus this ratio is extremely sensitive 
to the deviations of the nucleon DA from the asymptotic form. Therefore measurements of 
the neutron structure function in the near threshold region would considerably constrain the 
form of the nucleon distribution amplitude. 



Conclusions 



In this paper we derived a new soft pion theorem for the threshold pion production by a 
hard electromagnetic probe, i.e. with the probe of virtuality Q 2 3> A 2 (A ~ 1 GeV is a 
typical hadronic scale). This new hSPT allows us to express the pion production amplitudes 
in terms of the distribution amplitudes of the nucleon. The latter enter the description 
of various nucleon form factors at large momentum transfer. These new relations give a 
possibility to constrain further the nucleon distribution amplitude using data on threshold 
inelastic electron scattering from the nucleon at high momentum transfer. 

Here we restricted ourselves only to the case of the leading chiral contributions to the 
corresponding threshold amplitudes. The chiral corrections can be rather easily computed 
using methods of the chiral perturbation theory. 

Using a generic symmetric model for the nucleon DAs we demonstrated that various 
observables for near threshold pion production at high momentum transfer are sensitive 
to the parameters of nucleon DAs. This shows that the near threshold pion production 
by a hard electromagnetic probe is a new valuable source of information about nucleon 
distribution amplitudes. Studies with a broader range of models of nucleon DAs will be 
presented elsewhere. 
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Our analysis was restricted to the leading twist QCD contributions. The application of 
the methods developed here to the models for soft contributions to the baryon form factors 
(see a review in p(| [21]) would allow one to derive predictions of these models for hard near 
threshold pion production. This might be an exciting possibility to use hSPT to discriminate 
between soft and hard mechanisms for high momentum transfer reactions. 

We also note that the generalized ttN distribution amplitudes above the threshold are 
complex functions. Not very far from the threshold the phases of ttN DAs are fixed by 
known elastic ttN scattering phase shifts which allows us to use the Omnes solution to the 
corresponding dispersion relations in order to fix the dependence of ttN DAs on the mass of 
the corresponding nN system || . Since this dependence is fixed to a large extent by the 
ttN scattering phase shifts, the generalized ttN DAs also contain information about DAs of 
various nN resonances. In a complete analogy with the two-pion distribution amplitudes the 
generalized ttN DAs give a possibility to describe resonance and non-resonance contributions 
to the reaction (JJ) at high photon virtuality in a unified way. 

The formalism developed here can be also used for other hard reactions with the pion 
emission near the threshold. For instance one can analyze the Compton scattering 7 + N — > 
7 + (irN) at large momentum transfer. 

The study of the discussed processes should be feasible at the top of the current JLab 
energies and should be one of the high priorities of JLab at 12 GeV. 

We are grateful to P. Bosted, D. I. Diakonov, A. V. Efremov, L. Frankfurt, K. Goeke, 
V. Yu. Petrov, M. Praszalowicz, A. V. Radyushkin, N. G. Stefanis, P. Stoler and M. Van- 
derhaeghen for interesting discussions. This research has been supported by DFG, BMFB, 
A. v. Humboldt foundation and DOE. 
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